ROMEl/1423/06 
TIT/HEP-549 
|hep-th/0602123| 



Supersymmetry Breaking 

by Constant Boundary Superpotentials 
in Warped Space 



Nobuhito Maru^"^ ^, Norisuke Sakai*^'') ^, and Nobuhiro Uekusa^^^ 



(a) Dipartimento di Fisica, Universitd di Roma "La Sapienza" 
and INFN, Sezione di Roma, P.le Aldo Mow 2, 1-00185 Roma, Italy 

^''^ Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 



Abstract 

Supersymmetry breaking by constant (field independent) superpotentials local- 
ized at boundaries is studied in a super symmetric warped space model. We calculate 
the Kaluza-Klein mass spectrum of the hypermultiplet. We take into account of the 
radion and the compensator supermultiplets, as well as the bulk mass c for the 
hypermultiplet. The mass splitting is similar to that of the Scherk-Schwarz super- 
symmetry breaking (in flat space) for large |c|, and has an interesting dependence 
on the bulk mass parameter c. We show that the radius is stabilized by the presence 
of the constant boundary superpotentials. 
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1 Introduction 



Supersymmetry (SUSY) [T] is a well-motivated extension to the Standard Model, which 
plays a crucial role in solving the gauge hierarchy problem. Extra dimensions with flat 
space |21, [SI or with the warped space jlj are also an alternative solution to the gauge 
hierarchy problem. Considering both ingredients is natural in the context of the string 
theory and is often taken as the starting point in the phenomenological model of the brane 
world scenarios. In such a setup, we have to compactify extra dimensions and break SUSY 
to obtain realistic four-dimensional physics. One of the simple ways to realize it is the 
Scherk-Schwarz (SS) mechanism of SUSY breaking J^. 

It is known that the SS SUSY breaking is equivalent to the SUSY breaking by a (bulk) 
constant (field independent) superpotential in flat space [HIIZIIH]- These two scenarios 
generate the same mass spectrum. It is natural to ask whether this equivalence still 
holds in warped space. This issue has been discussed in several interesting papers so far 
jni CHI Ej ■ According to jH] , if one considers the gauging of a symmetry of the theory 
in warped space, which is completely broken by the boundary conditions, one is led to 
an inconsistent theory. This implies an inconsistency of the SS SUSY breaking in a 
SUSY Randall-Sundrum model. In (TU], they discuss the SS twist for SU{2)r in the five- 
dimensional gauged supergravity. Their conclusion is that if the background geometry is 
AdS4, SUSY is broken by the SS twist, but if the background geometry is the Randall- 
Sundrum geometry, SUSY is not broken by the SS twist. This statement agrees with that 
of [0]. Recently, [IJj also discussed this issue from the viewpoint of the five- dimensional 
conformal supergravity. It tells us that whether SUSY is broken by the SS twist or not in 
warped space depends on the way to gauge U{1)r. If we gauge U{1)r by the graviphoton 
with Z2-odd gauge coupling, SUSY is not broken by the SS twist. However, they pointed 
out that if we gauge f/(l)_R by the graviphoton with Z'2-even gauge coupling, it is possible 
to break SUSY. This seems to disagree with the statement in . In the light of these 
facts, the issue of the SS SUSY breaking in warped space is not settled. Furthermore, 
the equivalence of the SS SUSY breaking to the SUSY breaking by a constant (field 
independent) superpotential in warped space is also unclear compared to the fiat space 
case. A number of works have studied the SUSY Randall-Sundrum model fT^ ITHIIT^ IT^ . 

More recently, higher dimensional gauge theories are also used to provide Higgs fields 
as an extra-dimensional component of gauge fields jTHI. The true ground state of such a 
system is generally determined by the one-loop effective potential 5ZI-[2n]; which requires 
mass spectrum of all the Kaluza-Klein towers, especially their SUSY breaking pattern, 
since the effective potential vanishes in the supersymmetric limit. For phenomenological 
applications, it is also important to include bulk mass parameters for hypermultiplets 

m- 

The purpose of this paper is to investigate SUSY breaking effects and its properties 
caused by constant (field independent) superpotentials localized at fixed points in the 
SUSY Randall-Sundrum model. Taking the hypermultiplet and including the compen- 
sating multiplet and the radion multiplet consistently (from supergravity), we calculate 
the Kaluza-Klein mass spectrum and their mode functions of the hypermultiplet. We find 
that the mass spectrum depends on the bulk mass parameter in addition to the strength 
of the constant boundary superpotential. We observe a similarity of our result to the SS 
SUSY breaking spectrum in fiat space for large bulk mass parameter |c|. We also show 
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that the radius is stabihzed by the presence of the constant boundary superpotentials ^. 

This paper is organized as follows. The model is introduced in Sec|21 The mass 
spectrum of the hyper mult iplet is calculated in SecOland is compared with the SS SUSY 
breaking. The mode functions are obtained in SecEJ We show in SecISlthat the radius is 
stabihzed. SeclHl gives a conclusion. Appendices include a few details of calculations. 



2 Model 



We consider a five-dimensional supersymmetric model of a single hypermultiplet on the 
Randall- Sundrum background, whose metric is 



ds^ = e-^^'^rj^.^dx^dx'' + R^dy'^, a{y) = k\y\ 



(2.1) 



where R is the radius of of the orbifold S^/Z2, k is the AdS^ curvature scale, and 
the angle of is denoted by y{0 < y < n). In terms of superfields for four manifest 
supersymmetry, our Lagrangian reads (Hj 



£5 = J d'^e^if^if{T + T^)e~^^+^^^''{<^^<^ + <^'<l>^^ -6M, 



+ d'e 



$" 



a,- --C Ta' 



(2.2) 



where (y? = 1 + is the compensator chiral supermultiplet (of supergravity) , T = R + 
is the radion chiral supermultiplet, and $, are chiral supermultiplets representing 
the hypermultiplet. The Z2 parity is assigned to be even (odd) for $($'^). The derivative 
with respect to y is denoted by such as a' = da/dy. The five-dimensional Planck 
mass is denoted as M5. Here we consider a model with constant (field independent) 
superpotentials localized at the fixed points y = 0,n 



Wt = 2Mliwo5iy)+wJiy-7r)), 



(2.3) 



where Wo,Wt, are dimensionless constants which will be assumed to be C(l). 

Since interesting physics in the extra dimensions is contained solely in the part of the 
Lagrangian ()2.2|) containing auxiliary components F, F'^, Ft, F^, we extract that part 



-e-'^^" {2RF^ F + FtFV + i^F^^) 



2^'^(2i?0tF + Ft(0V - 3M|))(Ft - 4cr) + h.c.j + 
+e-2«'^i?(0V + - GMl){Fl - Fla){F^ - Fra) 



+ 

+e 



3e 

3 



ma 



{F^ - Ftct) 



-~c]Ra' 



F' 



--c\Ra' 



+ 



a, 



--c\Ra' 



c] Frcy'A + \i.c 



(2.4) 



*For related discussions on stabilization, see |21II22[0S| . for example. 
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where cj)^ and F'^ (0 and F) are scalar and auxiliary components of ($). The Lagrangian 
fl2.4|) gives the following equations of motion for auxiliary fields 





-Oy(p ' 


R 


1 


t J, Q J,ct 


' mf 








R 


1 










1 


R 


6M| 


3(1- 




•r 

r 






r 



2M| 



+ 



3M| 



6M| 



(2.5) 



2M| 



+ 



3M| 



(2.6) 



Mr 3Mr 6Mr ^ V2 / 2M^ 



3(1 -2i?(T)^^^ l-2i?a^^../3 
— -W}, H ^^^V^ ( - - c I i?cr' 



(2.7) 
(2.8) 



where the partial integration has been performed in (|2.5|) and r = 0^0 + (p'^'^ip'^ — 6M|. 
We can eliminate these auxiliary fields in ()2.4p by substituting ()2. 5^ - 1)2. 8|) . To obtain the 
mass spectrum, we take out only the part of the Lagrangian which is bilinear in scalar 
fields 



bilinear 
aux 



'r' 



-4Ra 



-d,<h'^+ - + c]Ra'<h'^ + 



2M3 



+ 



2M3 



+ 



m /3 
lM3 [2 



c Ra'chd)" + h.c. 



(2.9) 



The bilinear part of the Lagrangian turns out to be identical to that derived from another 
Lagrangian 



£5 = y"rf4^1y,t^(7- + 7-t)e-2iJ<x(^t^^^c^ct_g^^ 



dy-(--c]Ta' 



+ W\ + \i.c 



(2.10) 



which was proposed in Ref . [2^ based on a more accurate treatment of the radion superfield 
using super gravity. 



3 Mass spectrum of hyperscalar 

Let us calculate the mass spectrum of scalar component fields and 0"^ of the hyper- 
multiplet. To allow possible discontinuities of the Z2 odd field 0"^ across the fixed points 
?/ = 0, TT, we define 

0^(x, y) ^ e{y)h%x, y), e{y) ^ | _V<\''/o ' ^^"^^ 
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where h'^{x,y) is a parity even function with possibly nonvanishing value at y = 0,7r. 
The equations of motion for and (p'^ are given by 







2M| 



-2{S{y) - 6{y - n))h^ + + 7{e{y)fkRh^ - e{y)dyh^ 



2M| 



'e{y)h'^ 



(3.2) 



2M3 \^ 2M3 

2(a,%)-9,5(y-7r))/^^t + 



Wb + 29,0 1 - 2{26iy) - 25{y - n))dyh''^ 
dyWb 



2M3 



(3.3) 



One should note that the equations for and 0'^ couple only through the boundary 
superpotential Wb. In the limit of vanishing boundary superpotential, we should obtain 
a series of effective fields with 0(x, y) component only, and another series of effective 
fields with (j)'^{x,y) component only. For nonvanishing boundary superpotential, these 
two sets mix each other and we obtain n-th Kaluza-Klein effective field 0^(a;) with its 
mode functions 6^(?/) as 4>{x,y) component and b'^{y) as (j)'^{x,y) component 



0(X,?/) 

(f)''{x,y) 



/ n 1=1,2 \ 



Kiy) 



(3.4) 



where / is the indices corresponding to the two independent effective fields eigenvalues. 

Assuming that the effective four- dimensional field 0„(x) has mass m„, we easily find 
solutions in the bulk in terms of the Bessel functions J^, (J/3) and Fq,, (Yg) [TB] 



K{y) 
K{y) 



-,2Ra 



n 



J,(m„e^7A;) + 6,(m„)y„(m„e^7fc) 



Ra 





1, 


, a = 








, /5 = 


1 





(3.5) 
(3.6) 



where we have not yet specified which mode 7 = 1,2 until we determine mass eigenvalues 
later. 

Because of S{y),S{y — vr) in the boundary superpotential and also derivatives of sign 
function e{y), we obtain several types of singular contributions, S^{y) terms, dy6{y) terms 
and S{y) terms and similarly for y replaced hj y — n. To find out these singular terms, 
we use the following identity valid as a result of a properly regularized calculation ^ 



1 



S{ymy)y = 0^(2/), 5{y - nMy)y = -5iy - vr) 



(3.7) 



All these singular terms need to be canceled resulting in boundary conditions. The 
first boundary condition comes from 5^ terms in the equation of motion (|3.2p for 



- 26^(0) + wobniO) 
26^(7r) +w^6„(7r) 



0, 
0. 



(3.8) 
(3.9) 



''This follows from a wide range of regularization respecting the relation 26(i/) = de{y)/dy. 
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The second boundary condition comes from 6 function in the equation of motion for 

m„ 1 



7 1 
= -wobl{0)--wo 



- 46„(0) + 2g-c)6„(0)-^ 



7 X 1 

= - 3^7. 



Nr. 

mk'i 



{J'^inin/k) + baimn)Y^{mn/k)} 



(3.10) 
(3.11) 



+ 46„(7r) - 2 (| - c) 6„(7r) + ^ 



{ j;,(m„e^^7A;) + 6,(m„)r^(m„e«'=7fc)} 



with J'(-2) = dJ{z)/dz. We find that no additional boundary conditions ^ arises from the 
equation of motion of (f)'^, since the boundary condition from dy6 is identical to that from 

We wish to solve the boundary conditions (j3.8|) - (j3.1ip in the limit 



run/k <^ 1, 
mrfi^^^/k > 1. 



(3.12) 
(3.13) 



Because of fl3.12|l . the coefficient ha and hjj is small and Y^^Yp terms can be neglected 
at ?/ = TT. As given in Appendix |XJ the condition ()3.13p allows the remaining Bessel 
functions Ja and Jp to be approximated by their asymptotic forms. 
After changing variables, 



k 



2a + 1 
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-vr = a, 



2(3 + 1 



TT = b. 



(3.14) 



we can rewrite boundary conditions at y = vr ()A.10|) and ()A.11|) as summarized in Ap- 
pendix ^ into a matrix form 



" Mn M12 ■ 






M21 M22 







0. 



(3.15) 



Mn 
M12 

M21 
M22 



2 cos(x — 6), 
Wj, cos{x — a), 

— f (c^ — c + 1 1) cos(a; — b) + 2x sm(x — 6)) 
5 v / 

(c(l — c) cos(x — a) — 2x sin(x — a)) . 



(3.16) 
(3.17) 

(3.18) 

(3.19) 



Mass spectrum is determined by nontrivial solutions of ffuctuations require the van- 
ishing determinant of the matrix M 







2cos(x — b)[c{l — c) cos{x — a) — 2xsin(x — a)] 
9 , — c+11 , X 



—w^ cos a; — a) 



6 



cos(x ~ ^) + sin(a; — b) 



(3.20) 



'Other potentially singular contributions do not contribute because oi e x 5 — e 5^ — Q. 
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which determines the mass m„. It is most useful to start from the hmit of vanishing 
boundary superpotential, 10^,^0, giving two possibihties. The first solution is purely 
4>{x,y) modes {(p^i^x.y) = 0) 



= c(l — c) cos(a; — a) — 2x sin(x — a), 
which gives solutions for large x as 



X = niT + a + 



cil-c) 



O 



2 / ' 



2n-K 

with integer n. The second solution is purely (j)'^{x, y) modes (</'(x, y) = 0) 

= cos(x — b), 

which gives solutions 



X 



n - ^ ) TT + &, 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



with integer n. When ^ 0, these two modes ()3.22|) and ()3.24j) mix each other to form 
two independent modes with slightly different mass spectra and with associated mode 
functions. 

We can rewrite the eigenvalue equation (|3.2(jp in a more convenient form 



tan^(x — a) — A{x) tan(x — a) + B{x) = 0, 



(3.25) 



A{x)^^+(l + ^]cot{a-b), 5(x)^£cot(a-6)-^, (3.26) 



C = c{l-c) 



{c^ - c + 11)< 
12 



(3.27) 



Since we are interested in highly excited Kaluza-Klein states (n ^ 1), we can approximate 
x ~ nvr in the denominator of C/{2x). Thus we find 



tan(a; — a) 



A{n7c) 



1 ± 



\ 



ABimr) 



mx] 



(3.28) 



We can determine the sign choice in ()3.28|) by examining the solution near t/;^ ^ 0. Since 
the upper (lower) sign for < c < 1 (c < 0, 1 < c) reduces to the behavior in Eq. ()3.22|) at 
Wt, 0, it should correspond to the mode reducing to the purely 0(x, y) component at 
— *■ 0. We shall call this mode as the dominantly (p mode, and denote it by the mode 
functions (/ = 1) as bl^{y),b^{y) and the effective field as (plj"^{x) = 0„(a;) in equation 
()3.4|) . The other choice of the sign corresponds to the mode reducing to the purely (j)'^{x, y) 
component at — » , which behaves as in Eq. ()3.24|) . We shall call this mode as the 
dominantly (p^ mode, and denote the mode functions (/ = 2) as bl^i^y) , b'^ (y) and the 



effective field as 



X] 



[x) in equation (j3.4|) . 



Dominantly (p mode: 



a=i 
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We first consider the dominantly mode. Depending on the value of the bulk mass 
parameter c, there is a qualitative difference of mass spectrum as a function of w^. Let 
us first consider the case of |c| > 1/2. Then the parameter becomes a — b = ±7r/2, which 
makes the energy levels for purely and 0^ modes to be split only small amount of the 
order of 0{l/mT) at high enough excitation (n ^ 1). Therefore we need to distinguish two 
cases depending on the relative magnitude of \wt,\ and l/(n7r). For \wn\ ^ l/(?27r) ^ 1, 
we obtain perturbation from nondegenerate eigenvalues to find a deviation proportional 
to wl 



tan(x — a) 



c{l - c) 
2n'K 




(3.29) 



For ~ C'(l), we obtain perturbation from (approximately) degenerate eigenvalues to 
find a deviation linear in 



tan(x — a) ^ ± 
which gives the mass 



2v^ 



for l/2<c<l 

for c < -1/2 or c> 1 ' 



2a + 1 \ im 

n H : I TT ± 
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2v^ 



, ~ 0(1), 



(3.30) 



(3.31) 



where the plus (minus) sign should be taken for 1/2 < c< l(c< —1/2 orc> 1). 

We show tan(a; — a) with the upper sign (dominantly mode) as a function of Wt, in 

Fig.in 



tan(a; — a) 
0.4 




0.25 0.5 0.75 



1 .25 



Figure 1: The dependence of tan(x — a) with the upper sign (dominantly mode) in 
the case of |c| > 1/2. We use c = 0.5 and n = 10. 

For the case of |c| < 1/2, the parameter becomes a — b = cn, which makes the energy 
levels for purely and 0^ modes to be nondegenerate even at high excitation (n 3> 1). 
Therefore we always obtain a mass shift proportional to wl- For w-,, <^ l/(n7r), the 
solution is approximated by 

Wc(l — c) wl tan CTT 
x-a)^^ 2 . 3.32 
' 2mi 12 ^ ' 

For Wt, ~ C(l), it is approximated by 



,2\ / / ^2 
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This gives the mass whose w^^ dependence involves the bulk mass parameter c 

2a + 1' 



m„ ~ ke 



-Rkw 



< + 12 



n H ■ — I TT H ^ ' I 1 - W 1 + — tan^ cvr 

4 / 24tanc7r V 3 



(3.34) 



We show tan(a:; — a) as a function of w.,^ in Fig|21 

tan(a; — a) 



15 1 1.25 ™T 




Figure 2: The 1(77^ dependence of tan(x — a) with the upper sign (dominantly (p niode) in 
the case of |c| < 1/2. We use c = 0.1 and n = 10. 



Dominantly 0^ mode: 0^=^ = 

If we choose the lower sign in Eg. ()3.28j) . we obtain the dominantly 0^ mode. In the 
case of |c| > 1/2, we show tan(x — a) as a function of in FigOl For <^ l/(n7r) <^ 1 

ta,n{x — a) 




Figure 3: The w.,^ dependence of tan(x — a) with the lower sign (dominantly 0^^ mode) in 
the case of |c| > 1/2. We use c = 0.5 and n = 10. 



in the case of |c| > 1/2, we obtain 

tanfx — a 



niT 



6c(l-c) " 

For Wt, ~ C^(l) in the case of |c| > 1/2, it becomes 



tan(a; — a) ^ 



for l/2<c<l 



2\/3 ]+ foi' c<— l/2orc>l 



which gives the mass 



nir. 



ke 



-Rk-K 



' 1 , 2/5+ A \w^\ 

n TT =F — 1= 

. 2^ 4 i ^2v^ 



{n > 1) 



(3.35) 



(3.36) 



(3.37) 
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tan(;r- b + Ti/l) 




0.25 0.5 0.75 1 1.25 



Figure 4: The dependence of tan(x — 6 + 7r/2) for dominantly mode in the case of 
|c| < 1/2. We use c = 0.1 and n = 10. 

where the plus (minus) sign should be chosen for c< —1/2 orc> 1(1/2 < c< 1). 

In the case of |c| < 1/2, we show tan(a; — 6 + 7r/2) instead of tan(a; — a) in FigEJ It is 
approximated as 

, N wltarLcn 

tan(a; -6 + 7r/2) ^ ^ , (3.38) 

which gives the mass whose w-,^ dependence involves the bulk mass parameter c 

,2 



m„ ~ ke' 



-Rkn 



' 1 2/5 + 1\ <tanc7r 

n h — TT + — 

2 4/ 12 



(3.39) 



As for the mass spectrum of hyperfermions, one can immediately check that the lin- 
earized equations of motion for hyperfermions are not affected by the constant boundary 
superpotentials wq, Wt, locahzed at branes. Therefore the mass spectrum for hyperfermions 
is that for hyperscalars with tf^r = 0. 



4 Mode functions 

Let us determine the mode functions for each eigenvalue obtained above. The correspond- 
ing eigenvector of each solution gives the ratio Nn/N^ for the mode functions {pn^eb':^)^ 
in Eas. fpOll - dTT)!) . 

Dominantly (j) mode: 0^=^ = 

We first consider the dominantly modes corresponding to the choice of the upper 
sign in Eq. ()3.28|) . Using ()3.15|) - ()3.19|) . we find the eigenvector 

N- Mu M21' ^ ■ ' 

which is most easily evaluated by 

Mil = ^ \ , (4.2) 

Mil 2 cos(a — 6) — tan(a; — a) sin(a — h) ' 

For |c| < 1/2, we obtain the ratio Nn/N^ which is shown as a function of w-jr in Fig|31 
using the solution ()3.28|) with the upper sign for the dominantly modes. Using ()3.32|) 
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Figure 5: The ratio Nn/N^ as a function in the case of |c| < 1/2. We use c = 0.1 and 
n = 10. 



for small Wtt, it can be approximated as 



2 cos CTT 



(4.3) 



This result represents 0{wn) mixing of as we expected. 

For |c| > 1/2, we obtain the ratio by letting a — b = ±7r/2 in Eq. ()4.2|l . Using 
the solution ()3.28|) with the upper sign for the dominantly modes, we show Nn/N^ 
as a function of Wt, in FiglHl Using ()3.35|) for w <^ (l/nvr), the ratio Nn/N^ can be 



1.75 
1 . 5 
1.25 
1 

.75 
. 5 
0.25 



0.25 0.5 0.75 1 1.25 



Figure 6: The ratio Nn/N^ as a function 10.,^ for dominantly (f)'^ mode in the case of 
|c| > 1/2. We use c = 0.5 and n = 10. 



approximated by 







mrwT^ 


N^n 




c(l - c) 



(4.4) 



The result in this case represents 0{wn) mixing of b^ characteristic of the perturbation 
between nondegenerate states. Using ()3.36p for ~ C'(l), we obtain a wave function 
with the 0{1) mixing between and (f)'^ components 

§ « ^. (4-5) 

n 

as anticipated by the perturbation between the nearly degenerate {w-„ ^ l/nn) states. 
Dominantly mode: 0^"^ = (p'^ 
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Similarly we can obtain the dominantly (j)'^ modes corresponding to the choice of lower 
sign in Eq. ()3.28|) . It is now more convenient to use the ratio 



" - (4.6) 



22 



21 



(4.7) 



22 



(c^ - c + ll)tan(x - 6 + f) 



2x 



c(l — c)(tan(x — ^ + f ) cos(a — 6) — sin(a — h)) + 2x(cos(a — 6) + tan(x — 6 + |) sin(a — h)) 

For |c| < 1/2 case, we show N^/Nn as a function of in FiglTJ using the solution 
.28|1 with the lower sign for the dominantly (j)'^ modes. Using ()3.38jl for small 10-,^, the 




0.25 0.5 0.75 1 1.25 



Figure 7: The ratio N^/Nn as a function w-j^ for dominantly (f)'^ mode in the case of 
|c| < 1/2. We use c = 0.1 and n = 10. 



ratio is approximately given by 



n 



6 cos CTT 



(4.^ 



This result exhibits 0{wt^) mixing of (j) for the perturbation between nondegenerate 
states. 

For |c| > 1/2, we show N^/N^ as a function of in FiglHl using the solution (j3.28j) 
with the lower sign for the dominantly (j)'^ modes. Using ()3.35p for <ti I /(mi), the ratio 



-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 



0.25 0.5 0.75 1 1.25 



Figure 8: The ratio N^/Nn as a function w,r for dominantly 0^ mode in the case of 
|c| > 1/2. We use c = 0.5 and n = 10. 
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is approximately given by 









Nn 




3c(l - c) 



(4.9) 



which exhibits 0{wt,) mixing of as a perturbation of nondegenerate states. Using 
(|3.36p for Wj, ~ we obtain mixing between and 4>'^ components 



n 



1 

73' 



(4.10) 



as anticipated by the tu^ perturbation between the nearly degenerate {wt, ^ l/rni) states. 



5 Comment on the radius stabilization 

In this section, we study the background solution for hypermultiplets and possible radius 
stabilization in the presence of constant boundary superpotential Wb together with the 
compensator and the radion T supermultiplets. 

As a simple and realistic approximation, we consider a perturbative treatment for 
small values of boundary superpotential Wb- Let us first find out a background solution 
for the case of vanishing boundary superpotential: Wq = = 0. By requiring no SUSY 
breaking with vanishing auxiliary fields F = F'^ = = Ft = in Eqs. ()2.5j) - ()2.8|l . we 
obtain a unique supersymmetric solution with a complex constant N2, 



N2 exp 
= 0. 



-~c]Ra 



(5.1) 
(5.2) 



The potential vanishes for any values of the radius R and the complex parameter N2. 
Therefore the supersymmetric solution (pg has a flat directions along 0, as is typical for a 
SUSY solution. Moreover, radius R is undetermined in this Wq,Wt, limit. 

For nonzero wq and w-j^, we will find that supersymmetry is generically broken, and a 
nontrivial potential will be generated. Let us now assume |ifo| ~ I^ttI = w <S 1 and work 
out perturbative solutions of the equations of motion for x and x'^ as deviations from the 
supersymmetric solutions in Eqs. (j5.ip and (j5.2j) 



s + X, 
c 



ex 



(5.3) 
(5.4) 



where we allow possible discontinuities of the Z2 odd field 0^ across the fixed points 
y = 0,71, similarly to Eq. p.ip . The equations of motion from the Lagrangian ()2.4|) shows 
that 0^^^ arises in the first order in w 



+e 



- + c\Ra' 



-4Ra 



s'rs 



-ARa 



-6M| 



(?,-2cf{Ra'fex'- 



-dy + + c) Ra] 6X') 



'ARa 



2M| 



bl{dyWb 



0. 



(5.5) 



Provided the bulk fermion mass parameter c ^ 3/2, we can change a variable from y 
to a dimensionless variable X = (j)l(f)s/{QM^) — 1, and obtain the field equation for the 
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bulk {y ^ 0, vr) 



-{?>-2cYX\X + 1 



9X2 



-(3 - 2c)X(X + 1)(3 - 2c + 2(1 - 2c)X)^ 

+ (3 - 2cf + (3 - 2c) (9/2 - c)X + (1/2 - c)(3/2 + c)X'^ 



(5.6) 



whose solution is given for generic values of the bulk fermion mass parameter c (7^ 3/2, 1/2) 



as 



= X~\X + l)(5/2-c)/(3-2c) ^ ^^(^ ^ ^)-(l-2c)/(3-2c)(^ ^ (3 _ 2c)/(l _ 2c))] , (5.7) 

where Ci and C2 are constants of integration. The solution for c = 1/2 is given as 

= X~\X + 1) [c'l + c'2(X - ln(X + 1))] , (5.8) 

with another integration constants c'^^^c^- The field equation also yields the boundary 
conditions as the cancellation conditions of dy6{y) and dy6{y 



<Pl{dyW,) = 0, 



(5.9) 



respectively. Using Wb in Eq. ()2.3|) and the background solution (pg in Eq. ()5.ip . we obtain 
the boundary conditions more explicitly 



y=0 



y=n 



2X 



2X 



y=0 2(1 -iV)' 



y=TT 



2{Ne 



{3-2c)Rkn 



1) 



(5.10) 
(5.11) 



where we defined a dimensionless parameter N = |X2p/(6M|). These boundary condi- 
tions determine the integration constants Ci, C2 (or c[, c'2) in terms of the unique undeter- 
mined parameter N2. We find for generic values of c (7^ 1/2, 3/2) 

f \ ((l-2c)iVe2^'=^ + 2e-(i-2c)fifc^)^^^^l_2c)iV + 2)u;^e-^'=^ 
ci — — I — ^ 5-2c I ^ . — — ~ — ~ 5(5.12) 



^ 5-2e 
^2X2(3-2c) 



;i - 2c)X(e2-R'=- - 1) + 2(e-(i-2'=)-fJ*^- - 1) 



C2 



'1 - 2c) (wo + Wt^C 



,2iV#^ y (1 - 2c)X(e2«fc- - 1) + 2(e-(i-2c)R'=- - 1)' 
whereas we find the solution for c = 1/2 as 

, / iV| \ (iVe^^'^^ - In iV - 1 - 2Rkn)wo + (iV - In iV - l)ti;^e-^'=^ 



(5.13) 



,2X, 



]V(e2i?fc^ - 1) - 2Rk7r 



Wo + w^e 



-Rk-rr 



^2Nj N{e^R^^ - 1) - 2i?A;7r 



(5.14) 
(5.15) 
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In the case of c = 3/2, we can solve Eq. (|5.5|) in the bulk directly and find 



;^^ = c;'e^- + 4'e3^-, (5.16) 
whose integration constants determined by the boundary condition as 



2{N - 1) e2 
2(iV-l) e^^^^^ - 1 



From the equation of motion for 0, we find that the other perturbation x of 
background 0^ is of order 0{w'^), and that it is given in terms of the background 0^ and 
the first order perturbation ~ 0{w) uniquely, except for the amount of the admixture 
of the background solution 0^, which is always undetermined. Namely x is determined in 
terms of N2 without additional integration constants. 

By inserting these solutions into the Lagrangian (|2.4p and integrating over the extra 
dimension y, we obtain the potential as a function of the radius R and the complex 
normalization parameter N2 

y = ^ r dv( - 2c^iV^/2-2c+2/(3-2c)g((3-2c)(5/2-2c)+2)i?,CT 

2M| io I ^ 

+ + c + (3 - 2c) (-^ + 2c- [3(iVe(^-2^)^'^ - 1)]"')) ;^^t|^e-^^^ (5.19) 

where we should use the solution in Eq. ()5.7|1 subject to the boundary condition ()5.9|) for 
X'^. By performing integration and using the boundary conditions ()5.1()|1 and (j^.llj) . we 
find the potential 



V = -iVtA;c5iV5/2-2e+2/(3-2c) ^ ^^^ii3-2cf -2)Rkn^^ 

\N2\^kw?: / p 3-2c \ 

+ ^ — ^ -4c2 + 12c " 



4(1 - AT) V 3(1 -iV) / 

\l^h]jw^ _4^2 + _ 6 . ^ , (5.20) 

4(]Ve{3-2c)flfc. - 1) V 3(jVe(3-2c)R*- -1)/ 

where we should use C2 in Eq. (j5.13|) for generic values of c. It should be replaced by 
C2 in Eq. ()5.15|) for c = 1/2. We find that c = 1/2 case can be obtained as a smooth 
limit from c 7^ 1/2. We see that the supersymmetry breaking induced by the boundary 
superpotential produces the potential as a nontrivial function of N2 and R. 

Let us now study the stabilization of the radius R and the modulus For simplicity, 
we consider the case where = and the constant N2 is real. Namely we assume that 
the source of the SUSY breaking is localized on only one of the brane at y = (Planck 
brane). Then the potential becomes 



2 I (1 - 2c)(e2^'=- - 1)N + 2(e(2'=-i)«fe- 

N f o 3 - 2c \ 1 , , 

+ ^ -4c^ + 12c - 6 + ^ y 5.21 

1-A^V 3(1 -iV)/ J ^ ^ 
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We need to require the stationary condition for both modes R and N2 

dV dV 

— = and ^ = 0. (5.22) 
dR dN ^ ^ 

The former condition dV/dR = leads to 

N = e-(3-2'=)^fc-^ (5.23) 

whereas the latter condition gives 

2(4-2c- 3^)iV^~^^-^ 2(^e^Rkn _ i^ffi-2c-ji^^ 



jY(g2i?,fe7r _ 1) _ _A_(e(2c-l)Rfc7r _ l]Sl(^^2Rk7r - 1) - _2_ (e(2c-l)i?fe7r _ 



2 



+ f-4c^ + 12c - 6 + (^^Mi±ii) ) . (5.24) 

(i-iv)2v 3(i-iv) ; 

The vacuum expectation values for and {Rk) are obtained as solutions of Eqs. ()5.23p 
and (j5.24j) . depending only on the mass parameter c. We find that there is a unique 
nontrivial minimum with a finite value of the radius R and the normalization N2 for the 
flat direction provided c < Ccr with 



. (5.25) 

At the critical value of the mass parameter Ccr, the minimum occurs at infinite radius and 
vanishing normalization N2 

iV(cer) = 0, i?(c„) = 00. (5.26) 

To examine the stabilization for c < Ccr more closely, we parametrize c = Ccr — Ac with a 
small Ac. After using the relation Eq. ()5.23|l . N = g-(3-2c)RA;7r^ ^j-^^^^ potential 

(15.2111 for c = Ccr — Ac at the leading order of Ac and A^ consists of two pieces 

y . ?^(^. + ^,), (5.27) 

3 2 

F2 = -iV (^-8cer + y ) Ac. (5.29) 

The first piece Vi is positive and increases with a power larger than unity as a function of 
A^, whereas the second piece is negative and linear in N, whose coefficient is proportional 
to Ac. The potential V and its pieces Vi, V2 are depicted as a function of N in FigEl It is 
now obvious that a unique minimum occurs at finite values of N provided Ac > (c < Ccr) 
and that the minimum point approaches iV — > as Ac — (c — *• Ccr). Actually the FigEl 
demonstrates only the stability along the direction of N, after the other variable R is 
eliminated by the stationary condition ()5.23|) . We have checked that this minimum point 
gives a true minimum of the potential V{R, N) as a function of two variables, establishing 
the stability in both directions. For Ac 7^ 0, the stationary condition dV/dN = becomes 

^ ( 3 - 2ccr + 7^4^) N'-'"'''^^ - f ^ - 8cer) Ac. (5.30) 
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Figure 9: Potential for c = Ca — Ac. 



From this equation, we find that the stationary point at the leading order of Ac as 



R 



N 



[2(l-ccr)(3-2c„) + l] kTT 
(3 - 2c,,) (f - 4c, 



In 



(3 - 2c„) (f - 4c, 



l-ccr 



2(2cc 



1) 2 



l-ccr 
3-2ccr 



2(2Ccr 1) (2 Ccr 3-2ccr 

Ac 



Ac 



3-2ecr 
(3-2ccr){2-2ccr) + l 



Numerically R is evaluated as 



R 



10k 



In- 3.4 

Ac 



(5.31) 
(5.32) 



(5.33) 



which implies that Rk > 1 is satisfied for Ac < 10^®. Therefore our model with super- 
symmetry breaking by constant boundary superpotentials realizes the radius stabilization 
at Rk > 1 without requiring any additional mechanism. 

Another attempt for a possible approximation is given in Appendix O 



6 Conclusion 

In this paper, we have studied SUSY breaking by constant superpotentials at the bound- 
aries of the Randall- Sundrum geometry. We have calculated the Kaluza-Klein mass spec- 
trum of hypermultiplet in the SUSY Randall-Sundrum model including constant super- 
potentials localized at branes. In particular, we have taken account of the bulk mass 
parameter c, radion and compensator. 

Here we note that the presence of the compensator does contribute to the above 
calculated spectrum of hyperscalars. In the auxiliary field Lagrangian ()2.9|) . F component 
of the compensator gives the following type of terms 

Wb(l)dy(f)', Wb(f)'dy(j), W,a'<j)<f)', ly^^^V + ^^V^), (6.34) 

and their hermitian conjugate. From ()6.34|) . the equation of motion for includes 
Wf,dy(j)'^, W^(j)\ which are relevant to the boundary condition coming from the 5^ terms 



16 



as in (|3.9p and WbdyCp^, Wba'(f)^, which are relevant to the boundary condition coming 
from the 6 terms as in ()3.1ip . Therefore the boundary conditions ()3.15|) to determine 
mass eigenvalues is different if we do not introduce the compensator. 

In addition, we find that the SUSY breaking mass splitting and mode functions exhibit 
qualitatively different dependence for different values of the bulk mass parameter c. 

As mentioned in Introduction, it is not yet clear what is the SS SUSY breaking in 
warped space. However, it has been shown that the SS SUSY breaking and SUSY breaking 
by the brane locahzed constant superpotentials are equivalent in flat space [3 El IB]- 
Therefore our results should reduce to the SS SUSY breaking if we take the fiat limit. 
Since we have used the approximation m„/fc <^ 1 in order to focus on the effect of the 
warp factor, we cannot take the fiat limit immediately from our results. Nevertheless it 
is interesting to observe the following similarity with the SS SUSY breaking in our mass 
spectrum. Let us first note that the bulk mass parameter c should have large magnitude 
in order to take a proper fiat limit — > as seen from Eq. ()2.2|) . Our result for large |c| 
shows a linear dependence on w^, which is quite similar to that of the SS SUSY breaking 
as summarized in Appendix El This similarity seems to suggest the equivalence of the SS 
SUSY breaking with the SUSY breaking by constant superpotentials might be extendable 
in some way from fiat space to the warped space. This is in accord with a recent proposal 
that the SS mechanism works in the supergravity gauged by a Z2 even coupling ^Tj. On 
the other hand, we have found that the mass splitting depends on the bulk mass parameter 
for |c| < 1/2, which is a new pattern of SUSY breaking. It would be interesting that the 
different behavior for different values of the bulk mass parameter c might be understood 
by the analysis in the gauged supergravity. We hope that this work would give some 
insights to clarify whether the equivalence holds between SS breaking and SUSY breaking 
by constant superpotentials even in warped space. 

We have also discussed the possibility of radius stabilization. In the limit of vanish- 
ing constant boundary superpotentials, we have obtained a SUSY background solutions 
as a solution of the classical equations of motion. The solution has two fiat directions, 
the radius R and the complex moduli parameter N2 for the amplitude along one of the 
hypermultiplet scalar 0. Then, we have introduced constant boundary superpotentials 
wq^tt <^ 1 as small perturbations. For the simplest situation of = (boundary super- 
potential only for the Planck brane), we obtain the deviations from the SUSY solutions 
in the leading order of wq. Using these solutions, we have explicitly calculated the radion 
potential, and have shown that the radius R and the moduli N2 are indeed stabilized at 
finite values. 

When the classical solution is trivial (j) = (j)'^ = 0, the radion potential is generated at 
1-loop after SUSY is broken and the radius might be stabilized by the Casimir energy. If 
SUSY is broken at high energy, the radion is stabilized for a special bulk mass parameter 
c = 1/2 25 . This mechanism might be helpful in the present case. 

The radius stabilization has been studied also in the AdS4 background where SS SUSY 
breaking can be formulated. In models with nonzero superpotential it has been 
found that hypermultiplets give positive contributions to the radion potential, contrary 
to the negative contributions from the gravity multiplet. This provides various patterns 
of radion potential. However explicit calculation with arbitrary values of c remains to 
be studied. It would be interesting to apply our analysis for arbitrary values of c to the 
radius stabilization of the AdS4 background. 
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A Solving boundary conditions for mass spectrum 

The boundary conditions from 5"^ terms ()3.8p and (j3.9|) can be rewritten in terms of the 
Bessel functions 



[Ja{mn/k) + ba{mn)Ya{mn/k)] 



[J/sirrin/k) + hp{mn)Yp{mn/k)] , (A.l) 



A„ 



2 



w^A„= 



J^(m„e^'^7A;) + 6^(m0y^(m„e^'=7A;) . (A.2) 



Similarly, the boundary conditions from 5 terms ()3.10p and ()3.11|) can be rewritten 

'5 
—'i 
6 



Wq-1- 2c j (Jairrin/k) + 6a(m„)Y'„(m„/fc)) 

iJ'ai^n/k) + ha{mn)Y^{mn/k)) + \wj^^ \j'Jnin/k) + hi3{nin)Y'Jnin/k) 

3 K I\„ 



k 



(A.3) 



- 1 - 2c) (j,(m„e^'=VA;) + 6„(m„)V;(m„e^'=VA;)) 



2m. 



n ^Rk-n 



1 rUn Nn 



k 

--W. 



3 ' k Nf^ 

We can determine ba, hp from the boundary conditions (jA.lll and ()A.3j) at ?/ = 

1 



(A.4) 



{p + 5)wi - 6(1 + 2c - 2a) 

(-(/5 + 5K + 6(l + 2c + 2a) 
1 



X 



Jaimn/k) ^^^^ NnJpimn/k) 



Ya{mn/k) 



N^Y^{mn/k) 



(/? + 5)u;g - 6(1 + 2c - 2a) 

+ 5K-6(l + 2c + 2a)) (^) 



X 



2a 



TT 



\2kJ r(a + l)r(a) 



-4/?wo 



vr 



N^\2kJ r(/? + i)r(a) 
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(A.5) 



hp{mn) 



(/3 + 5)u7§ - 6(1 + 2c - 2a) 



X 



12awo 



J^(m„/fc) at;; 

Yp{mn/k) Nn 

1 



+ ((/?-5V^ + 6(l + 2c-2a)) 



Jp{mn/k) 
Yp{mn/k) 



(/3 + 5)u7g-6(l + 2c-2a) 



X 



-12«^o ( ^ 



a+/3 



77 



V{a + l)r(/5) iV„ 



((/?-5H^ + 6(l + 2c-2a)) (^)'' 



TT 



2A;; r(/? + i)r(/5) 



(A.6) 



where we used the approximation for |2;| <^ 1 



2, 



r(/3 + i) 



- I sin(7ra)r(l — a) 



TT 



■ (A.7) 



In the approximation (|3.12p . we can neglect since {mn/ky°',{mn/k)°'~^^,{mn/k)'^^ 

are very small. Using the other approximation condition ()3.13|) . we can approximate the 
Bessel functions for \z\ » 1 as 



Ya{z) 




2 

nz 
2 

TIZ 



2a + I 

cos I Z ■ TT 




2a + I 

Sm I Z ■ TT 



Aa'^-l 
8z 



sm z — 



2a + 1 



-TT 



2a + 1 
sm I z ■ — vr 



(A.8) 
(A.9) 



Therefore, it is a good approximation to impose the boundary conditions at y = tt with 
only Ja{i3) terms disregarding Ya(/3) terms. Then, the approximated boundary conditions 
at y = TT (|3.1(J|) and (|3.11|) can be read 








-— cos( — — e 
k 



2/3 + 1 nin RkTT 2a + 1 



(A.IO) 



N, 

+ 



,171, 



c(l — c) cos( — — e 
k 



n ^Rkw 



2a + 1 



2m. 



4 



n RkT, 



,m. 



, -e smi — — e 

k K 



2a + 1 



-vr 



n 



& -C+ll 



cos( — — e 
k 



n „RkiT 



2(3 + 1 



it) + 



Rk-K 



,m 



n Rkn 



2/? + l 



vr) 



(A.ll) 



which gives the matrix eigenvalue equation (|3.15p . 



B Relation to SS SUSY breaking 



In this appendix, we briefly review that the SS SUSY breaking and SUSY breaking by 
the brane locahzed constant superpotentials are equivalent in fiat space O 13 E] and 
compare our spectrum in the SUSY Randall-Sundrum model with that in fiat space. Let 
us consider the following hypermultiplet Lagrangian in fiat space 



1$ 



c\2\ 



j d^e^^dy^ + h.c. 



(B.l) 
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(B.4) 



Its auxiliary field Lagrangian can be read as 

= (l^r + \Ff) + - F<9,0= + h.c] , (B.2) 

then, it is easy to solve 

Ft = dy(f)', F^^ = -dy(f). (B.3) 
Here, we twist the hyperscalar in terms of R-symmetry rotation SU{2)ii, 

(^t )(l/ + 2-i?) = e---(^t )(^)- 
where q is SU{2)ii charge. This can be written as 

(;,) fe) |, ) (!/). (B.6) 
where 0^ are periodic functions for and f{y) was found in [7], 

f{y) = — I — ^y) + — I — viy) (b.6) 

where 

ri{y) = 21 + 1, l7cR<y < {1 + 1)'kR, {I : integer) (B.7) 

is the "staircase" function that jumps by two units for every ttR along y. Note f{y + 
271 R) = f{y) + Wq + Wt^, which implies Wq + Wt, = 27rg to yield a correct twist. Then the 
equation of motion for (f) is 

{-R'v'^d.d. - - _ f'^dy4>'^ _ f~^) = 0. (B.8) 
The equation of motion for 0^ is given by 

{-RY'd.d, - dl)4>' + dy{f4>^) + f'{dy4>^ + r~^'^) = 0. (B.9) 
These equations can be solved by 

4) = cos{ny + f) (B.IO) 

4)^ = sm{ny + f) (B.ll) 

where f{y) is the periodic function given by 

r m^m±^(^2l7T-y) for 2/7r < y < (2/ + l)7r 

j^J v^^m±^(^(^2l + l)n -y) for (2/ + l)7r < y < 2(/ + l)7r (B.12) 
[ for y = Itx 

f = w,5{y) + wj{y - vr) - !f^L±^. (B.13) 
The mass spectrum is obtained as'' 



1 

m„ = — 
R 



n . (B.14) 

27r 

This spectrum in flat case is similar to our results and fl3.37j) in the Randall- 

Sundrum model for large |c| since only the scalar fields in the hypermultiplet receive the 
mass shift which is linear in w^. 



^Here, a regularization different from ours in H3.7|l is used for S{y)e'^{y) by taking the limit for e{y) 
first without respecting the relation 2d{y) = de{y)/dy. This treatment results in S{y)e^{y) = 6{y), instead 
of our 5{y)e^{y) = 5(y)/3 as in (PTTjl . 
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C Another approximation for background solution 



In this appendix, we consider another approximation to explore solutions of the equations 
of motion. 

We see that F components are nonlinear in t^*^, resulting in a difficulty to solve 
the equations of motion for 0, 0*^. To examine the leading order effects of the constant 
boundary superpotential, we attempt the following approximation 

0,</)^ ~ 0{R~'/^), wo,.r^O{l), (C.l) 
R-^ < Ms, i?(T ~ C(107r). (C.2) 

Evaluating the order of magnitude of the terms in the equations of motion for auxiliary 
fields fl2.5j) - ()2.8|l . the dominant terms are found to be 



-Ra 



F 



R 

-Bo 



2M| 



R 



2M| 



(C.3) 
(C.4) 



R 



1 



-Ra 



ml'"""'"' 



6Rcr ^4 a - 1 3 ( 1 + i?or) - 2Rac 



6M| 



(C.5) 



2 / MV 



+ 



Ml 



(C.6) 



Then we find that the auxiliary fields Lagrangian is exactly given by the bilinear La- 
grangian (j2.9p . Let us look for the classical solution for the background because the 
radion potential can be obtained from the classical solution. We assume the classical 
solution for the background to depend only on y. Since we obtain the bilinear Lagrangian 
()2.9|) . the classical equations of motion is identical to the linearized one in ()3.2j) and ()3.3|) . 
except that dfj^d^ term in the last lines should be discarded. 
We can easily solve the equations of motion in the bulk 



,(3/2-c)i?,<7 
,(3/2+c)iJ(T 



[^Cg(l-2c)i?,a ^ 



2J, 
2], 



(C.7) 
(C.8) 



where Ni^2,N^2 integration constants. The boundary conditions from 5^ terms (j3.8|) 
and (j3.9|) become 



= -2{N'^ + N^)+wo{Ni + N2), 

The boundary conditions from 6 terms ()3.10|) and (j3.1H) become 

"1 /3 \ _ 1 

.3 



(C.9) 
(C.IO) 







-Wo 



g + c) (iVr + iV2^) + (l-2c)^iVr 



7 



+ -Wo(iVr + iV2)-2iVi(l + 2c), 



-(- + c) (Ar{=e(i-2^)^'=" + N^) + (1 - 2c)-N^e^^-^'^'^''^ 
3 V 2 / 3 



+ -(iV{=e(^-2'=)-f^'=" + N^) + 2iVi(l + 2c)e 
3 



Rkn 



(C.ll) 



(C.12) 
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From (fTTTTIl and (frn2|) . we obtain 



N2W0 



Rkn 



X 



12(1 + 2c) + 9/2 + c 

(w, + u;oe^^'") /(9 + 2c) 
- (w^e(i-2c)ijfc^ ^ w;oe^'=^) /(II - 2c) 



-\ -1 



11/2 



(C.13) 



From 



9|) and ()C.13|1 . we obtain 



N2W0 



+ 



^^g(i-2c)flfc^ + Woe 



Rkw 



12(1 + 2c) 9/2 + c 

xw7oW^,(l - e(i-2c)m.)/(^2(l + 2c)). 



11/2 -c 



(C.14) 



It is interesting to see that iVf , A^|, Ni ^ for — ^ 0. We observe that the overall 
constant N2 remains undetermined, because of the bilinear Lagrangian in our approxima- 
tion. Moreover, we do not have good reasons to assume this constant N2 not to depend 
on the radius R, since all these normalization constants A^^i, A^2, N^, N2 are related by the 
boundary condition involving R. Therefore we cannot give a definite answer on the radius 
stabilization in the present approximation. 

If we take into account higher order corrections and/or additional mechanisms, we 
can presumably determine the normalization constant, although it is a very hard task to 
analyze. 
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